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Quantum invariants of the orbit dependent pairing problem are identified in the limit where the 
orbits become degenerate. These quantum invariants are simultaneously diagonalized with the help 
of the Bethe ansatz method and a symmetry in their spectra relating the eigenvalues corresponding 
to different number of pairs is discussed. These quantum invariants are analogous to the well known 
rational Gaudin magnet Hamiltonians which play the same role in the reduced pairing case (i.e., 
^s^j , orbit independent pairing with non degenerate energy levels). It is pointed out that although the 

^ ' reduced pairing and the degenerate cases are opposite of each other, the Bethe ansatz diagonalization 

of the invariant operators in both cases are based on the same algebraic structure described by the 
rational Gaudin algebra. 



Oh- 

C3 



o 



Keywords: Nuclear Pairing, Bethe Ansatz, Exact Solution, Quasispin, Gaudin Algebra, Integrability. 



I. INTRODUCTION 



Strong pair correlations are observed in fermionic many body systems which energetically favor large wave 
function overlaps. This phenomenon is known as pairing and it plays an important role in our understanding 
of many body physics (See Ref. [l| for a review). Historically, the physical significance of pairing was 
first realized with the microscopic theory of superconductivity developed by J. Bardeen, L. N. Cooper and 
J. R. Schrieffer (BCS) in 1957 [2]. Following the success of the BCS theory, the idea of pairing was carried 
over to other areas of physics as well. In particular, pairing now plays an essential role in the nuclear shell 
model as the residual interaction between nucleons and successfully recounts for various properties of atomic 

oo ' nuclei 

In order to investigate the influence of pairing on nuclear properties, many authors have used exact 
analytical solutions of nuclear shell model which are available in some simplified cases. For example, in Ref. 

Kerman considered pairs of nucleons coupled to angular momentum zero occupying a single orbit and 
introduced the quasi-spin formalism in which these pairs can be treated within suitable representations of 
the angular momentum algebra. He used this formalism to write down the exact energy eigenstates and 
to analyze the influence of pairing on the collective vibrations of nuclei. Quasi-spin formalism can also 
be extended to the case of several orbits in which case the quasi-spin angular momenta corresponding to 
different orbits commute with each other and the pairing term has the form of a coupling between these 
angular momenta. This observation establishes a direct link between the fcrmion pairing models and the 
interacting spin models (see Refs. g 0, | for reviews). The exact diagonalization of the later model 
was carried out by R. W. Richardson in 1962 in the case of the orbit-independent (i.e., reduced) pairing 
interaction Q. Later, it was clarified that the exact solvability of the pairing Hamiltonian in the reduced 
pairing case can be understood in terms of a set of quantum invariants which commute with one another and 
also with the Hamiltonian 
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111 ]. These invariants are called rational Gaudin magnet Hamiltonians since 
they stem from the work of M. Gaudin who was originally trying to find the largest set mutually commuting 
operators for a given system of interacting spins. In Ref. 10], Gaudin also showed that the pair creation 
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and annihilation operators which are used in building the simultaneous eigenstates of the rational Gaudin 
magnet Hamiltonians form an algebra which is today known as the rational Gaudin algebra. It is worth 
mentioning that the rational Gaudin algebra is related to the rational solution of the classical Yang-Baxter 
equation which appears as an integrability condition in many contexts. As a result, the rational Gaudin 
magnet operators and the rational Gaudin algebra have found many other applications in physics (see Refs. 
@, [J 3 for reviews and Refs. Q, [T3I . EH , Ea . Ell for some interesting applications). They have also been 



generalized to include other underlying algebraic structures (i.e., higher rank algebras, super-algebras and 
deformed algebras) besides the angular momentum algebra. There is an extensive literature on this subject 
and the interested reader may find the Refs. 0, El, El El |23| | useful. 



Although the Richardson-Gaudin solution is successfully used in nuclear physics, the assumption of reduced 
pairing sometimes proves to be too stringent. In many cases, the effective residual interactions between the 
nucleons are best described by a pairing force whose strength differs between the orbits 1 . The Hamiltonian 

^ = EE £ J fl J ra a i m ~ \ G \J2 C i° 3> a ] m a ] - m a j> m' C 1 ) 

is frequently used to describe such an orbit dependent pairing interaction. Here, j denotes the total angular 
momentum of an orbit and ej denotes its energy. The overall strength of the pairing term against the kinetic 
term is measured by the constant \G\ which has the dimension of energy whereas the relative pairing strengths 
are measured by the dimensionlcss constants Cj. Richardson-Gaudin scheme mentioned above applies to the 
special case of this Hamiltonian in which all Cj's are the same whereas all single particle energy levels Ej 
are different from one another (the reduced pairing case). The focus of this paper, however, is the opposite 
case in which all the Cj 's are different from one another and all single particle energies Ej are the same (the 
degenerate case). 

The problem described by Hamiltonian given in Eq. |T]) is exactly solvable in both the reduced pairing 
case and the degenerate case. As mentioned above, in the reduced pairing case the solution was given 
by Richardson-Gaudin scheme and the corresponding quantum invariants are the rational Gaudin magnet 
Hamiltonians. In the degenerate case, the exact energy eigenvalues and eigenstates were obtained in a series 



of papers by Pan et al 25| and by Balantekin et al 26l.l27| and the purpose of the present paper is to identify 



the corresponding quantum invariants in the degenerate case. In addition, it will be shown that the quantum 
invariants in the degenerate case can be simultaneously diagonalized with the help of the algebraic Bethe 
ansatz method. An interesting observation regarding the Bethe ansatz diagonalization is that although the 
reduced pairing and the degenerate cases are opposite of each other, the Bethe ansatz diagonalization of the 
invariant operators in both cases is connected with the rational Gaudin algebra. 

The organization of this paper is as follows: Section II is a brief review of the quasispin formalism and it 
also serves to introduce some notation. In Section III, a short review of the Richardson-Gaudin formalism 
and the rational Gaudin magnet Hamiltonians is presented. The main results of this paper, i.e., the quantum 
invariants in the degenerate case and their simultaneous diagonalization with the Bethe ansatz method are 
presented in Section IV. This section also contains a discussion about a symmetry in the spectra of these 
quantum invariants relating eigenvalues corresponding to different number of particles. In Section V, we 
consider the rational Gaudin algebra and point out its relationship with the Bethe ansatz diagonalization 
in both the reduced pairing and the degenerate cases. Section VI summarizes the main conclusions. The 
details of some of the Bethe ansatz calculations can be found in the Appendix. 



1 See, for example, the treatment of the neutron pairs occupying the valance shell of Ni isotopes by Auerbach 
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II. QUASI-SPIN FORMALISM AND THE EXACT SOLUTIONS OF THE PAIRING 

HAMILTONIAN 



In the quasi-spin formalism, nucleoli pairs coupled to angular momentum zero are created and annihilated 
at the level j by the operators 

S+=E (-I)'" 7 ""] rrA -m «T = £ (-1)^*% - m aj m , (2) 

m>0 m>0 

respectively. Together with the operator 

Sj = — ^aj m a.j m + Clj _ m dj -m — (3) 

m>0 

they obey the well known angular momentum commutations relations 

[S+, Sp] = 25 jf S? [Si Sp] = ±5 3r Sf. (4) 

As a result, one has an angular momentum algebra (the so called quasi-spin algebra) for each orbit j such 
that those angular momenta corresponding to different orbits commute with one another. The pair number 
operator for the orbit j is given by 



Nj = - ^2 U] m a 3 m + a] _ m a,j _ m ) . (5) 



2 

m>0 



It is related to the operator 5° given in Eq. © by the formula 



S° = N 3 -^ (6) 



where ilj is the maximum number of pairs which can occupy the level j. Note that j is always an half integer 
because of the spin-orbit coupling in the nuclei. As a result, if there are no unpaired particles at the level j, 
then 

^ =3 + \- (7) 

Also note that the pairing term in the Hamiltonian given in Eq. |T]) does not act on the unpaired particles. 
If there is an unpaired particle at the level j, its effect will be i) to add a constant Sj to the Hamiltonian 
because of the kinetic term and ii) to reduce the maximum number of pairs which can occupy the level j 
by one, i.e., to take Qj to Clj — 1. But here it will be assumed that there are no unpaired particles in the 
system. In this case, Eq. (jSJ) implies that 

_ ^ < s o < 2, (8) 

i.e., quasi-spin algebra corresponding to the level j is realized in the fij/2 representation. Therefore, in 
addition to the physical angular momentum quantum number j, we also have the quasi-spin quantum number 
£lj/2 for each level. The states 

|^,-^) and |^,^) (9) 

respectively represent the situations in which i) the level j is not occupied by any pairs and ii) it is maximally 
occupied by pairs. In the presence of several orbits with angular momenta j\, . . . , j n , the state 

10} = |%, 6) |%, 6) • - ® |%, (10) 
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represents a shell which contains no pairs whereas the state 

10} = 1%, %> ® 1%, %) ® • • • ® |%, %) (11) 
represents a shell which is fully occupied by pairs. 

The pairing Hamiltonian given in Eq. |T]) can be written in terms of the quasi-spin operators given in Eqs. 
© and © as 

£ = 5> (2S° + tlj) \G\(£ C ^+)(E c ^)- ( 12 ) 

3 3 j' 

Note that the operator ^ . e,-Sj in Hamiltonian (fT2|l creates a pair of particles in such a way that Cj can be 
viewed as the probability amplitude that this pair is found at the level j. For this reason the coefficients Cj 
are usually called occupation probability amplitudes and they are normalized as 

Es 2 = i- (13) 

3 

Although an occupation probability amplitude is a complex number in general, the parameters Cj can be 
taken as real without loss of generality. Because if one Cj is complex, a unitary transformation can always 
be performed on the quasi-spin algebra corresponding to the level j to make that Cj real. Also note that the 
Hamiltonian in Eq. (|12p contains a constant term 2ejQj which comes from Eq. This constant term 
is not dropped because it guarantees that the energy of the empty shell is zero, i.e., 

H\0}=0. (14) 

Using the commutators given in Eq. I©, one can show that the fully occupied shell |0) is also an eigenstate 
of the Hamiltonian with the energy 

H\0)=J2 (2^ - |G| Cj 2 ) fij |0>. (15) 

3 

Unlike the empty shell |0) and the fully occupied shell |0), the eigenstates of the pairing Hamiltonian 
corresponding to a partially occupied valance shell are unknown in the most general case. But, as mentioned 
in the Introduction, exact energies and eigenstates are known in the two opposite cases. Namely, the reduced 
pairing case characterized by 

Sl ^ S2 ^ • • • ^ e n 

ci = c 2 = • • • = c„, (16) 

and the degenerate case characterized by 

£l = £2 = ■ ■ ■ = £n 

ci ^ c 2 ^ ■ ■ ■ ^ c n . (17) 

These solutions will be reviewed in the next two sections together with the corresponding quantum invariants. 
But before closing this section, mention must be made of a third case in which exact eigenstates of the pairing 
Hamiltonian are known. This solution is available in the presence of two orbits with unequal energies and 
unequal occupation probability amplitudes, i.e., 



£1 ^ £2 

ci + c 2 . (18) 

But this third case will not be considered in this paper. Because the main interest of this paper is the 
quantum invariants of the pairing Hamiltonian and in the case of a two level system we have only two 
quantum invariants which are simply the Hamiltonian itself and the total pair number operator. 
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III. REDUCED PAIRING AND THE GAUDIN MAGNET OPERATORS 



In the reduced pairing case, described by Eq. (|16p . the pairing Hamiltonian given in Eq. (|12p becomes 

H R =J2^(^ + ^)-\G\dJ2s+S-. (19) 

3 3,3' 

Here d = 1/n is known as the level spacing and its appearance in the Hamiltonian is due to the normalization 
condition (fT3|) . Using a variational technique, Richardson showed in Ref. [9( that the eigenstates of the 
Hamiltonian given in Eq. (I19|) containing N pairs of particles are in the form 



J + (Zi)J + -(6)... J + (Zn)\0) (20) 

where the pair creation operators J + (£) are given by 

j 

and |0) is the state with no pairs defined in Eq. (JTUJ) . The values of the parameters £i, £2, ■ ■ ■ , £,n which 
appear in Eq. (f2T))) are to be determined by solving the system of equations 

^^l^-^+^O 1 ^ (22) 

simultaneously for k = 1, 2, . . . , N (see Ref. 9}). These equations generally have several distinct solutions. 
For each one of these solutions we have an eigenstate in the form of Eq. ([20]) and the corresponding energy 
is given by 

N 



£W=X> (23) 
fe=i 

The quantum invariants of the Hamiltonian in Eq. (|19[) are the rational Gaudin magnet Hamiltonians 



mentioned in the introduction 



10. 



ll| . They are given by 2 



R 3 =S°-\G\dY: 4^ ( 24 ) 



£7" — 

i'&j) 3 3 
where S j ■ S f is defined as 

• ~g s , = S*S° + ~ (S+S- + S-S+) . (25) 

The rational Gaudin magnet Hamiltonians mutually commute with one another and with the Hamiltonian 
H R , i.e., 

[R j ,R j ,]=Q [4,^=0 (26) 

for every — 1,2, ... , n. The Hamiltonian itself is not an independent invariant and it can be written in 
terms of the operators Rj as 

H R = J2 ~ \°\ d ) k i + \G\dY,RjRr - \G\dJ2^s ■ ~Sj + E £ i<V ( 27 ) 
j 3 j 



2 Strictly speaking, the operators studied by Gaudin himself did not include the one body term S^. 



G 



Similarly, the total pair number operator can also be written in terms of the operators Rj as 

* = E (** + t) (28) 

3 

As a result of Eq. , the eigenstates of the Hamiltonian are at the same time simultaneous eigenstates 
of rational Gaudin magnet operators as well. Let us denote the eigenvalues of the invariant Rj corresponding 
to the eigenstate with N pairs by Ej . In other words, 

4|0)=^ 0) |0) (29) 

for the empty shell and 

Ri J + (a) ■ • • J + {Zn)\Q) = eW J+(ti)J + fa) . . . J + (6v)|0) (30) 

for the eigenstates containing N pairs described by Eqs. flUgg) . The eigenvalues Ef ] and E^ N) are given 

by 

(o) % \G\d \^ fyfy 

jW) 3 3 



and 



£?w=Ef +iG|dx;2- s v (32) 



respectively. 



The pairing Hamiltonian in Eq. (|19p is only one of the exactly solvable models which can be built 
using the rational Gaudin magnet Hamiltonians. Various other linear or nonlinear combinations of rational 
Gaudin magnet operators can be used to built other useful exactly solvable models (see, for example, Refs. 

113 EJ). 



IV. INTEGRABILITY IN THE DEGENERATE CASE 



In the case of several orbits having the same energy but different occupation probability amplitudes, i.e., 
when the conditions in Eq. (|17|) are satisfied, the first term in the pairing Hamiltonian given in Eq. (|12|) 
becomes a constant which is proportional to the total number of pairs in the shell. Discarding this term, 
one can write the Hamiltonian as 

Hn = -\G\C£c j Sf)C£c f Sr l ). (33) 

3 j' 



Exact eigenvalues and eigenstates of the Hamiltonian given in Eq. ([55]) were obtained in Refs. [25l |26|. l27j. 
The purpose of this paper is to introduce the corresponding quantum invariants, i.e., the set of operators 
which commute with one another and with the Hamiltonian in Eq, 



The fact that the rational Gaudin magnet operators given in Eqs. (|24p mutually commute with one another 
is independent of the values of the parameters Sj . Naturally, one can try to replace the parameters Sj in the 
Gaudin operators with some arbitrary functions of Cj and try to determine the form of these functions so 
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that the new operators commute with the Hamiltonian in Eq. (|33|) as well. It turns out, however, that such 
as course of action does not yield the quantum invariants of the Hamiltonian given in Eq. (133p 3 . 

In order to find the invariant operators one can consider general number conserving operators in the form 
/' = AjS] + B 3 S+SJ + D rr Spl + F n' ( S t S ? + S Pi) ( 34 ) 

where Aj, Bj Djj> and Fjj> are some arbitrary coefficients. The condition that the above operators commute 
with one another and with the Hamiltonian in Eq. (|3"5|) gives us the allowed values of these coefficients. A 
straightforward calculation shows that the desired operators are given by 



Pj - S+Sj + 2 ^ 2 J _„2 S j S j' + E 73^772 ( S f S j> + S f S j 



(35) 



These operators mutually commute with one another 



P Pi 



= 



(36) 



for every j and j'. They also commute with the Hamiltonian given in Eq. (|33|) and with the total pair 
number operator: 



Pj,H D 







P h N 



0. 



(37) 



The Hamiltonian Hjj and the total number operator TV are not independent invariants but they are related 
to the operators Pj by the formulas 



\G\Y, c2 A=Hd 



(38) 



and 



Y t p i = N 2 -NC£ t ^ i + i) + \ E n J°j'- 



(39) 



3,3 



As a result of Eqs. (|3^|) and l|37p . the invariants Pj have the same eigenstates as the pairing Hamiltonian 



Ho given_in Eq. (|33|) . These eigenstates were given in Refs. 
method 



25 



2G 



27j with the help of the Bethe ansatz 



29| . In what follows, the corresponding eigenvalues of the invariant operators Pj will be presented. 



A summary of the results of this Section can be found in Table |U 
Following Refs. 



25 



2G 



27j . let us introduce the pair creation and annihilation operators 



-cf-x 



q+ 

2„ D 3 



S~(x) 



q- 



(40) 



Here, x is a complex variable and are the quasispin operators introduced in Eq. |J5J|. The Hamiltonian 
in Eq. ([55]) itself can be written in terms of these operators as 



-|G|5 + (0)5 _ (0). 



(41) 



3 In addition to the rational Gaudin magnet Hamiltonians, Gaudin also studied the so called trigonometric and hyperbolic 
Gaudin magnet Hamiltonians which also mutually commute with one another. Similar to the rational case, however, one 
cannot choose the parameters of trigonometric or hyperbolic Gaudin magnet Hamiltonians so as to make them commute with 
the Hamiltonian in Eq. Il33jl . 



The eigenstates of the Hamiltonian Hjj which are also simultaneous eigenstates of the invariants Pj can 
be written in terms of the pair creation and annihilation operators in Eq. (|40j) . Below, these eigenstates 
which were obtained in Refs. 25j, 1261 . |2jj will be reviewed in the order of increasing number of pairs and the 
corresponding eigenvalues of the invariant operators Pj will be given. 



Empty shell: The empty shell |0) given in Eq. (|T0|) obeys 



P J |0) = K (o) |0) 



(42) 



where Ej is given by 



E 



(0) 



(43) 



Eigenstates with N = 1: The eigenstates with one pair of particles fall in two classes. The state 

S+(0)|0) (44) 

is an eigenstate where S + (0) is obtained by putting x = in the operator given in Eq. (|40f . This state was 
first suggested by Talmi in Ref. and was shown to be an eigenstate of a class of Hamiltonians including 



the Hamiltonian in Eq. 

In addition to the state in Eq. (|44|) . the state 

S+(x)\0) (45) 
is also an eigenstate if a; is a solution of the Bethe ansatz equation 

y ^ML = o. (46 ) 

^ l (?. -x v ' 

j J 

The eigenvalues of the operators Pj corresponding to the eigenstates described above will be denoted by 
A^ and (Xj , respectively, i.e., 

PjS+mO) = \Vs + (0)\0) (47) 
PjS+ixM = S+(x)\0). (48) 

The eigenvalues A^- 1 and fj,^ can easily be computed using the commutators given in Eq. ([4]) together with 
Eqs. ([40]) and (j46| as follows (see the Appendix) 

Af) = Ef - P,, (49) 
^ = (50) 

J J 1 — CjX 

Note that the eigenstate in Eq. (|44p is unique whereas the eigenstate in Eq. (|4"5)) represents several eigen- 
states. Because in general the Bethe ansatz equation (|4l)]) has more than one solutions and for each one 
of them we have an eigenstate in the form of Eq. (|45|) . As a result, Eq. (|48|) also represents several 
eigenvalue-eigenstate equations . 

Eigenstates for 2 < N < N max /2: The results given above can be generalized to the states corresponding 
to a shell which is at most half full. Let N max = f^- denote the maximum number of pairs which can 
occupy the shell in consideration. Then for 2 < N < N max /2 the results obtained for one pair generalizes 
as follows: The state 

S + (0)S + (z 1 )...S + (z N . 1 )\0) (51) 
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which has N pairs of particles is an cigenstate if the parameters Zk are all different from one another and 
obey the following system of Bethe ansatz equations 

-%/2 1 ^ 1 



=-+ E > ( 52 ) 

J- / Cj Z m Z m k—1 ^ m ^ 



for every m = 1, 2, . . .N — 1. In addition, the state 

S + ( Xl )S + (x 2 )...S + (x N )\0) (53) 

which also has N pairs of particles is an eigenstate if the parameters Xk are all different from one another 
and satisfy the following system of Bethe ansatz equations: 

~"/ /2 = J2 - , (54) 

3 k=l 
(fe^m) 



for every m = 1,2, ...,N. Note that the states given in Eqs. ([51]) and (J53J) can be thought of as the 
generalizations of the states in Eqs. (|44|) and (|45)) . respectively. 

The Bethe ansatz equations given in Eqs. and ([54")) have more than one solutions. 4 Each one of these 
solutions gives us an eigenstate in the form of Eqs. (fSTj) and (|53|) . respectively. One should keep in mind, 
however, that the states in Eqs. (151j) and (|53j) are invariant under the permutations of the parameters and so 
are the corresponding Bethe ansatz equations. Consequently, solutions of the Bethe ansatz equations which 
differ only by a reordering of the variables should be counted as one solution. 

The eigenvalues of the operators Pj corresponding to the eigenstates given in Eqs. ([BTj) and (|53|) above 
will be denoted by Xj and [A , respectively, i.e., 

P 3 S + (0)S+( Z1 ) . . . S+(z N ^)\0) = Af > S+(0)S+( Z1 ) . . . S+(z N ^)\0) (55) 
Pj S+(x 1 )S+(x 2 ) . . . S+(x N )\0) = S+( Xl )S+(x 2 ) . . . S+(x N )\0) (56) 

These eigenvalues can be computed with the help of Eqs. (|4]), (|40j) . (f52|) and ([54]) . The results are given by 
(see the Appendix for the details) 

N-l „ 



N 



„f > = Ef> - T (58) 



Eigenstates for N max /2 < N: In order to write down the eigenstates and eigenvalues corresponding to 
a shell which is more than half full, we introduce the operator 



4 It is worth emphasizing that in driving Eqs. H52I I. all the parameters zj. were assumed to be different from each another and 
similarly for the Eqs. H54I I. As a result, any solution which contradicts this assumption should be discarded. In fact, it can 
be shown that if any two parameters in the state l|51|l are equal to each other, then this state cannot be an eigenstate. A 
similar result is also valid for the state in Eq. (153 D - 



TABLE I: Summary of the energy eigenvalues and the eigenstates of the pairing Hamiltonian in the degenerate limit, given in Eq. (|33[1 and its quantum invariants given 
in Eq. ((35} . Here, N max denotes the maximum number of pairs which can occupy the shell. 
Empty shell: N = 



State 



Bethe Ansatz Equation 



Eigenvalue of Pj 



Eigenvalue of Hd 



|0> 



No BAE 



E (o) _ Sb. 



~2~ y 1 ■■■■■ 



c/c, 



One pair of particles in the shell: N = 1 



State 

S + (x) |0) 
S + (0) |0> 



Bethe Ansatz Equation 



^ l/c? -X = 



I/. 

No BAE 



Eigenvalue of Pj 



E 



(0) 



1 - <* x 



Eigenvalue of Hd 



At most half full Shell: N < N max /2 



State 



S + {x 1 )S + (x 2 )...S + (x N )\0) 
S + (0) S + (zi)...S+ (z N -i)\0) 



Bethe Ansatz Equation 



. l/c? ~x k + ^ 



= 



Eigenvalue of Pj 

AT 



E 



i=i 



x k - XI 



+ E 



1/c? - z k z k " z k — Zl 
J 3 (=1 
(Mfc) 



(0) 



fc=l 2 
JV-1 

°i - E I 



Eigenvalue oi Hd 



\G\J2cPi + \G\J2 



Zk 



More than half full Shell: N max /2 < N 



State 



S~ (zi)S- (z 2 )...S~ (z N -i)\0) 



Bethe Ansatz Equation 



^7 2 v — ^ 2 _ 

/"a" 2 + — +> =° 

/c . - Zfc z fe f— f 2fc — z ; 



Eigenvalue of Pj 

JV-l 



(0) 



Eigenvalue of Hd 



\G\J2cPi + \G\J2 



Zk 



Full Shell N = N n 



State 
10) 



Bethe Ansatz Equation 



No BAE 



Eigenvalue of Pj 



E^ - a< 



Eigenvalue of Hd 
-|G|E C ^ 
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following Ref. 27]. This operator transforms the empty shell given in Eq. (fTU)) into the fully occupied shell 
given in Eq. (|11|) . i.e., 

ft|0) = |0). (60) 

In addition, it also transforms the pair creation operators into pair annihilation operators and visa versa: 

f^Sff = Sf f *<S°f - -S° (61) 

Therefore, the operator T* transforms a state containing N particle pairs into states containing N hole-pairs. 
However, there is no particle-hole symmetry in the problem described by the invariants Pj as can be easily 
verified by showing that 

Pj,f]?0. (62) 
On the other hand, it is easy to show that the operator 

B = f^S~(0) (63) 

commutes with the invariants Pj : 

[P j7 B]=0 [H D ,B]=Q. (64) 

This tells us that if \ip) is a simultaneous eigenstate of the operators Pj so is B\ip) unless \ip) is annihilated 
by B. It is easy to see from the definition of B that if the state \ip) has N particle-pairs, then the state B\tp) 
has N — 1 hole-pairs. Because the operator B first annihilates one pair and then replaces the remaining 
N — 1 particle-pairs with hole-pairs. Consider, for example, the eigenstates with one pair of particles given 
in Eqs. flH| and fl5]). One can easily show, using Eqs. I[4T]). pi)) and flU]) that 

BS+(0)|0) oc |0) (65) 

where |0) is the state which is maximally occupied by pairs defined in Eq. (jlip . Consequently, the states 
5 I+ (0)|0) and |0) have the same eigenvalues given by Eqs. (|4"5|) and (|50[) . On the other hand, the state in 
Eq. (|4"5|) which also has one pair of particles is annihilated by the operator B: 

BS + (x)\0)=0 (66) 

This can be easily verified using the commutators (UJ) and the Bethe ansatz equation (14"1)|) . Similarly, if we 
act on the eigenstate with N pairs of particles (2 < N < N max /2) given in Eq. (|51|) . we find that 

B S + (0)S + ( Zl ) . . . S + (z N ^)\0) ex S~(zi) . . . 5-(^-i)|0) (67) 

As a result the state 

^-(2ri)...^-(^_i)|0> (68) 

which has N— 1 pairs of holes, has the same eigenvalues as the state in Eq. (|51l) . i.e., those given in Eqs. ([57]) 
and |58p. On the other hand, the state in Eq. (|S"3")) which also has N pairs of particles (2 < N < N max /2) 
is annihilated by the operator B 

B S + { Xl )S + {x 2 )...S + (x N )\0) =0 (69) 
as can be verified by using the commutators ^ and the Bcthc ansatz equations given in Eqs. (|54p . 



The eigenstates and the eigenvalues of the invariants Pj described in this section are summarized in Table 
HI This table also contains the corresponding eigenvalues of the pairing Hamiltonian Hp given in Eq. (f3"3"]) 
obtained from the eigenvalues of its invariants by using the formula (|38p . For example, the eigenvalue of Hd 
corresponding to the state in Eq. (|4"4"1) is given by 

ICIEST (70) 

3 

etc . . . The energy eigenvalues of the Hamiltonian Hp obtained from Eq. (|38p are in agreement with those 

'life 



obtained earlier in Refs. 



25. 
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V. RATIONAL GAUDIN ALGEBRA 



Rational Gaudin algebra naturally appears from the pair creation and annihilation operators used in 



building the simultaneous eigenstates of the rational Gaudin magnet Hamiltonians [10L 1 1 71 ] - Formally, it 
is defined as an infinite dimensional algebra whose generators J + (A), J~ (A), J°(X) depend on a complex 
variable A and obey the commutation relations 



[J+(A),J-( M )] = 2- 



A — // 

[J°(X),JHH)} = ± J± ^- J± ^\ (71) 

A — /I 

for A 7^ fi. Commutators of J + (A), J - (A) and J (A) at the same value of the complex parameter are given 
by taking the limit fx — > A. We are usually interested, however, with the finite dimensional (unfaithful) 
realizations of this algebra. For example, it is not difficult to show that for any set of real parameters otj 
which are all different from one another, the operators 

J+K;A) = £^ J^X)^^ ^;A)=£^ (72) 

3 3 3 

form a realization of the rational Gaudin algebra given in Eqs. (|71[) . Note that in Eqs. ([72]) we explicitly 
stated the dependence of the generators on the parameters aj . For different values of these parameters, one 
obtains different realizations of the rational Gaudin algebra. 



The realization of the rational Gaudin algebra given in Eqs. (|72|) is related to the pairing problem in 
both the reduced pairing case and the degenerate case for different values of the parameters aj . By taking 
aj = 2ej in Eq. ([72"]) . one obtains the operators given in Eq. (|2"Tj) which are used in building the eigenstates 
in the reduced pairing case, i.e., 

J+(2 £i ;A) = £^-^ j- (2£j . ; A)=£^-_ jO (2£ .. A) = £_jl_ (73) 



The pair creation operators given in Eq. (|40|) can also be obtained from a realization of the rational 
Gaudin algebra although a change of basis is necessary. By taking aj = 1/cj in Eqs. ([72"]) one obtains 

g+ go 
J + Uc3-A) = Y,y^X J -^) = J2jJ^- X Al/^;A) = Ei7^A (?4) 

3 3 3 

and then the operators in Eq. (|40j) can be written as a linear combination of these operators, i.e., 

s±{x) = JHl/^yffl + JHl/cf-VZ) (75) 

Rational Gaudin algebra defined in Eqs. ()71|) appears in connection with various integrable models as 
mentioned in the Introduction. It is often the case that the step operators which are used in building the 
eigenstates live in one of its realizations. Here we see that the integrability of the pairing Hamiltonian in the 
reduced pairing and the degenerate pairing cases can be studied in two different realizations of the rational 
Gaudin algebra, namely those given in Eqs. ([73]) and ([74]) . respectively. 



VI. CONCLUSIONS 



In this paper, we have obtained the quantum invariants of the pairing Hamiltonian in the degenerate case 
and simultaneously diagonalized them with the help of the algebraic Bethe ansatz method. Although in 



13 



order to obtain the eigenvalues one should first solve a system of Bethe ansatz equations which are nonlinear 
and coupled to each other, solving them usually proves to be much more convenient then a direct numerical 
diagonalization method. Exact analytical methods for solving the Bethe ansatz equations also exist in some 



simplified cases (see, for example, Refs. lg, l26l |31|). The quantum invariants obtained in this paper are 
the counterparts of the well known rational Gaudin magnet Hamiltonians which play the same role in the 
reduced pairing case. It is worth mentioning that since the quantum invariants are mutually commuting 
operators, they can be used to build various other integrable models besides the ones considered in this 
paper. 

We also pointed out that the integrability of the pairing Hamiltonian in both the reduced pairing and the 
degenerate cases is connected with the rational Gaudin algebra. The generalizations of the rational Gaudin 
algebra to different underlying algebraic systems (such as higher order Lie algebras, quantum algebras or 
superalgebras) have been used to study the reduced pairing model and the related rational Gaudin magnet 
Hamiltonians in more general frameworks. The question then naturally arises whether or not one can do the 
same generalization for the degenerate pairing model and the related quantum invariants too. For example 
can the invariant operators of the degenerate pairing given in Eq. (|33|l be generalized to other algebraic 
systems and then used to study different integrable many body systems? The answer of this questions goes 
beyond the scope of this paper and will be considered elsewhere. 



APPENDIX A: OBTAINING THE EIGENSTATES WITH BETHE ANSATZ METHOD 



The simultaneous eigenstates and eigenvalues of the degenerate the operators given in Eq. (j3"5|) can be 
obtained using the method of algebraic Bethe ansatz. In this method, one first constructs a Bethe ansatz 
state 29] which includes some undetermined parameters and then substitutes this state into the eigenvalue- 
eigenstate equation Pj\ip) = Ejip). The requirement that the Bethe ansatz state obeys the eigenvalue- 
eigenstate equation yields a set of equations called the equations of Bethe ansatz, whose solutions determine 
the values of the parameters in the Bethe ansatz state. For example, in order to obtain the eigenstates with 
one pair of particles, one can start from a generic state in the form 

5+(x)|0) (Al) 

where S + (x) is defined in Eq. (|4()p. Using the commutators given in Eq. Q, one can show that the action 
of the operator Pj on such a state is given by 

P^+(z)|0) = (V - -%-) S+(x)\0) + [ ] ^2-|0). (A2) 



Clearly, if we choose x in such a way that the second term on the right hand side of Eq. (| A2|) vanishes, i.e., 
if 

4^ = o, (A3) 



/ Cj 2 - X 



then Eq. (|A2[) becomes an eigenvalue-eigenstate equation. One way is to satisfy Eq. (|A3|) is to take x = 
in which case Eq. (|A2j> yields the eigenvalue-eigenstate equation given in Eq. (jlT)) . Alternatively one can 
choose x in such a way that it satisfies the Bethe ansatz equation 



Qj = (A4) 



L/c^ -X 

in which case, Eq. (|A2ft yields the eigenvalue-eigenstate equation given in Eq. ((48 
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These results can be easily generalized to a state in the form 

S + (xx)S + (x 2 )...S + (x N )\0) (A5) 

which has N pairs of particles where N < N max /2. The parameters Xi,x 2 , . . . ,xn are in general complex 
and they are all different from one another. In fact, by acting on the state given in Eq. (|A5|) with the 
operators Pj one can easily show that if any two of the two parameters x±, x 2 , . . . , Xjy are the same, then 
the state in Eq. (|A5|) cannot be an eigenstate. If, on the other hand, the parameters Xi,x 2 , . . . , xn are all 
different from one another, then by acting on the state given in Eq. (IA5jl with the operators Pj given in Eq. 
(|35|. we find 

P.S+ (a*) S+ (x 2 ) ... S+ (x N ) |0) = f Ef - n % x \ S + (a*) S+ (x 2 ) ... S+ (x N ) |0) (A6) 



k=l 



+ X\ 



+ x 2 




+ 



+ XN 



N 

E 



, , x N -x k <—'l/c 2 . l -x N 

k=l y J 



S+ (x 2 )...S+ (x N )\0) 



■ S+ (x N )\0) 



S + ( X1 )S + (x 2 ).. ° jS J 10). 

1 — CjXpf 



Clearly, Eq. (|A6[) becomes an eigenvalue-eigenstate equation for all the operators Pj if we choose the variables 
xi,x 2i . . . ,xn so as to satisfy 



■T-k 



N 

E 



0,4 







(A7) 



for every k = 1, 2, . . . , N. 

If the variables x\,x 2 ,. . . ,xn are all different from zero, then they clearly have to obey 

2 ^ fL 



N 



= 



(A8) 



Xk-xi '-f i/cj - x k 
and in this case, Eq. (IA6|) yields the eigenvalue-eigenstate equation given in Eq. |56|) . 

If, on the other hand, one of the parameters, say x\ , is chosen to be zero (we cannot choose more than one 
Xk to be zero since the parameters must be different from one another) then the Bethe ansatz equation (|A7|) 
is automatically satisfied for fc = 1. In this case, the remaining parameters x 2 , . . . ,Xn are to be found by 
solving the N — 1 equations 

N 



E 



l/Cj — Xfc £fc * — ' Xk ~ Xi 

J 3 l=2{l^k) 







(A9) 



for A: = 2, 3, . . . , AT. If we rename the remaining N — 1 variables as X2 = Zi, = Z2, ■ ■ ■ , xn = zjv-i, then 
Eq. (|A6|) becomes the eigenvalue-eigenstate equation given in Eq. (f55]> . 
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